Quantum Smoluchowski equation for driven systems 



OS 

o 
o 

(N 
< 



I 

C 

o 
o 



> 

o 
in 
o> 
o 

o 

o> 
o 



X 



Raoul Dillenschneider and Eric Lutz 

Department of Physics, University of Augsburg, D-86135 Augsburg, Germany 

(Dated: April 21, 2009) 

We consider a driven quantum harmonic oscillator strongly coupled to a heat bath. Starting 
from the exact quantum Langevin equation, we use a Green's function approach to determine 
the corresponding semiclassical equation for the Wigner phase space distribution. In the limit of 
high friction, we apply Brinkman's method to derive the quantum Smoluchowski equation for the 
probability distribution in position space. We further determine the range of validity of the equation 
and discuss the special case of a Brownian parametric oscillator. 
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Introduction. The investigation of Brownian motion 
in the overdamped regime is of fundamental importance 
both from a theoretical and experimental point of view. 
In the limit of high friction, the velocity part of the phase 
space distribution of a Brownian particle quickly relaxes 
to equilibrium. Then for times much larger than the ve- 
locity relaxation time, a description in terms of the posi- 
tion distribution alone becomes possible and the Kramers 
equation reduces to the simpler Smoluchowski equation 

[Ell- 

Recently, a quantum generalization of the Smolu- 
chowski equation for systems strongly coupled to a heat 
bath has been put forward [3, 4]. The quantum Smolu- 
chowski equation provides a semiclassical description of 
the evolution of the diagonal matrix elements of the sys- 
tem density operator in the coordinate representation; it 
allows to study the influence of thermal as well as quan- 
tum fluctuations. The quantum Smoluchowski equation 
has been applied to problems involving both undriven 
and driven systems; examples include the study of the 
quantum decay rates for driven potential barriers [5j, 
quantum phase diffusion and charging effects in Joseph- 
son junctions jfj, quantum diffusion in tilted periodic po- 
tentials 0, [1| , and quantum extensions of nonequilibrium 
fluctuation theorems [9j. However, the exact expression 
of the quantum Smoluchowski equation for undriven sys- 
tems is still the subject of discussions 11, 12, 13, 14 1 
and a rigorous derivation for driven quantum systems is 
lacking. 

The aim of this paper is twofold: We first present 
a transparent derivation of the quantum Smoluchowski 
equation that makes use of controlled approximations, in 
order to complement the mathematically involved path 
integral derivation of Ref. Q and the heuristic approach 
of Ref. [12] . Second, we provide the first derivation of the 
driven quantum Smoluchowski equation and determine 
its range of validity. In the following we consider a quan- 
tum harmonic oscillator with arbitrary time-dependent 
frequency. We employ a simple Green's function ap- 
proach to determine the semiclassical expression of the 
diffusion coefficients appearing in the evolution equation 
of the Wigner phase space distribution, starting from 
the known exact quantum Langevin equation. In the 



large friction limit, we then apply the method introduced 
by Brinkman [2?| to derive the quantum Smoluchowski 
equation. In the undriven case, we recover the equa- 
tion obtained in Refs. 12 1 and [14J. In the driven case, 
on the other hand, we show that the quantum Smolu- 
chowski equation remains valid provided the driving rate 
is smaller than the velocity relaxation rate. We finally 
apply our results to the important example of the Brow- 
nian parametric quantum harmonic oscillator (l5| . 

Quantum Langevin equation. The starting point of 
our discussion is the quantum Langevin equation for 
a harmonic oscillator with time-dependent potential, 
V(x,t) = muj 2 (t)x 2 /2, linearly cou pled to a bath of in- 
dependent harmonic oscillators ID, 17 1, 



mx - 



dt 7(* 



t) X {t) + - 



-j(t)x + F(t). (1) 



In the above equation the damping kernel is given by 
j(t — t ) = 1/7T dv [J(v)/v) cos vit — t ), where J{v) 
is the spectral density of the bath. In the sequel we focus 
on the Ohmic regime where J(v) = m^v, the parameter 
7 denoting the friction coefficient. The fluctuating force 
operator F(t) verifies the correlation function 

(F(t)F(t)+F(t')F(t)) = 

°° dv , , (I3hv\ 

— J(i/)coth cosi/(t-t)> ( 2 ) 

-oo V 2 J 

where the average is taken over the initial bath degrees 
of freedom. Equation ^ is exact and is derived under 
the assumption that the initial density operator of sys- 
tem plus bath factorizes. The solution of the Langevin 
equation {1} can be written as, 



x(t) = mx Gx(t) + mx Q G 2 (t) + X(t), 



(3) 



where Xq and Xq are the initial position and velocity of 
the quantum system and X(t) = J Q dt G(t, t )F(t ) is the 
fluctuating position operator. The two functions G\{t) 
and G2(t) are solutions of the homogeneous Langevin 
equation ([T]) with vanishing member on the right hand 
side. By introducing the two solutions, </>i(£) and 4>2(t), 



2 



of the equation y + VL 2 (t)y = with Vl 2 (t) = uj 2 {t) -7 2 /4, 
we can explicitly write G\{t) = (1/m) exp(— 74/2) <f>i{t) 
and G 2 (4) = (l/m)exp(-74/2)<^ 2 (4). The Green's func- 
tion G(4, t ) is then given by the combination of the two 
functions and <p 2 taken at different times t and t , 

G (t,t) = ^ (fc (*)&(*') - <h(t')Mtj) ■ ( 4 ) 

TO V / 

A detailed derivation of the Green's function (H|) for 
the time-dependent harmonic oscillator can be found in 
Ref. [Hj]. The two functions 4>i(t) and <^> 2 (4) are linearly 
independent and their Wronskian obeys </>i(4)</) 2 (4) — 
(f>\(t)<j)2(i) = 1 for all time t. By fixing the initial con- 
ditions, <f>i(0) = <fa(0) = and 4> 2 (0) = <Pi(0) = 1, the 
Green's function G(4, 4 ) verifies the relations G(t, t) = 
and G (1 ' 0) (4,4) = -G^ ' 1 ^*,*) = 1/m, where the num- 
bers in parenthesis denote derivatives with respect to the 
first (second) time argument. Moreover, since the sys- 
tem is initially decoupled from the heat bath, we have 
(x Q F(t)) = (x F(t)) = 0. 

Wigner function. The phase space dynamics of a quan- 
tum system is conveniently described using the Wigner 
quasiprobability distribution [l9| . The evolution equa- 
tion of the Wigner function W(q,p, t) for a harmonic os- 
cillator coupled to a bath of harmonic oscillators is of the 
general form [U l^l , 



d 

—W 

dt 



■^W + m Q 2 (t)q-^-W + 2T(t)^-( P W) 
771 oq Op op 



-D„{t)^W + D„(t)^W. 



(5) 



An expression for the time-dependent parameters fi 2 (4), 
2T(t), D pp (t) and D qp {t) for a driven harmonic oscilla- 
tor can be directly derived from the quantum Langevin 
equation |T]) by generalizing the method introduced in 
Ref. (23)]. The first step is to rewrite Eq. (p} in a form 
that is local in time [24[ • This can be done by first invert- 
ing Eq. ([3]) and its first time derivative in order to express 
the initial coordinates xq and xq in terms of the time- 
dependent variables. The latter coordinates are then in- 
jected into the second time derivative of Eq. (|3|), leading 
to a Langevin equation with time-dependent coefficients, 



x + 2T(t)x + n 2 (t) = F(t)/m. 



(6) 



Here we have defined the function 2F(4) = 
(Gi(4)G 2 (4) - G 2 {t)G 1 {t))/L{t) and the parameter 
n 2 (t) = (G 2 (t)G 1 (t) - Gi(t)G 2 [t))/L(t). The denomi- 
nator is given by L(t) = m(G 2 (4)G x (4) - Gi(4)G 2 (4)). 
In the Ohmic regime, we simply have 2F(4) = 7 and 
£t 2 (t) = u> 2 {t). The second step is to derive an equa- 
tion similar to Eq. ([6]) from the evolution equation 
of the Wigner function (JSJ) by evaluating the first 
moments of the position and momentum operators, 
(x) = IZ.^IZoMq + (ih/2)d p )W(q,p,t) and 



(p) = m(x) = f^dqj^dpfr - (ih/2)d q )W(q,p,t). 
This yields 



(x.) + 2T(t)(x) +n 2 (t)(x) = 0. 



(7) 



Since the Langevin equation and the equation for the 
Wigner function describe the same process, the average 
of Eq. ((6]) must be equal to Eq. ([7]), implying that the 
time-dependent parameters 2T(t) and f2 2 (4) are identical 
in both equations. By repeating the same argument for 
the second moments of the position and momentum op- 
erators, the diffusion coefficients D pp (t) and D qp (t) can 
be related to the fluctuatingposition operator X(t) and 
the noise operator F(t) via 22, 23|, 



D PP (t) = -(X(t)F(t)+F(t)X(t)), 
D qP (t) = l(X(t)F(t)+F(t)X(t)). 



(8) 



Semiclassical diffusion coefficients. The diffusion co- 
efficients D pp (t) and D qp (t) can be further expressed in 
terms of the correlation function ([2]) of the noise opera- 
tor and the Green's function ([4} by using the definition 
of the fluctuating position operator X(t). We find, 

D PP (t) = £ f dtG^\t,i){F{t)F{t)+F{t)F{t)), 



D qP {t) = - J dt G{t,t){F{t )F(t) + F(t)F(t )). (9) 

It is important to notice that the Green's function G(4, t ) 
is the same for quantum and classical oscillators and that 
the quantum-mechanical nature of the process is solely 
encoded in the noise correlation function. A semiclas- 
sical approximation of the diffusion coefficients can ac- 
cordingly be obtained by expanding the noise correlator 
((2|) in powers of K. From Eq. ([2]) , we obtain up to second 
order for an Ohmic bath, 



(F(t )F(t)+F(t)F(t )) 



h2 imf3 



s'\t 



5(4 - 4 ) 

t') + o(h 3 ) 



(10) 



The first term on the right hand side corresponds to the 
classical noise correlation function, while the second term 
accounts for the first quantum corrections. We next eval- 
uate the classical and quantum contributions to the diffu- 
sion coefficients D pp {t) and D qp (t), Eqs. ([9]), separately. 
The classical expression of the coefficient D pp (t) reads 



&'&W(t,t')^8{t-t) = ??l, (11) 



6 



where we have used the relation G' 1,0 )(4, t) — 1/m. Since 
G(t,t) = the classical coefficient D qp (t) vanishes, 

D c qp {t) = \f dt'G(t, t')^S(t-t) = 0. (12) 
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In the limit h — > the diffusion coefficients are thus con- 
stant and Eq. ([5]) reduces to the familiar classical Klein- 
Kramers equation 0, Q . 

The quantum expressions of the diffusion coefficients 
on the other hand are obtained by considering the quan- 
tum corrections to the noise correlation function (fTU)) . 
For the coefficient D pp (t), we have, 



D q pp (i) 



m 

T 



dt 



jh 2 m(3 



27m 2 A ( ., 



*(*-*), 

(13) 



where we have used G^ 1,2 \t,t) — —(us 2 (t) — 7 2 )/m and 
introduced the parameter A = H 2 f3 2 /24m [12|; the latter 
measures the magnitude of quantum fluctuations. By 
noting moreover that G^°' 2 '(t,t) = — 7/m, we find that 
the quantum contribution to D qp (t) is 



DUt) = - 



dt 



'G{t,t){- 



■yh 2 m(3 n f.\ 2j 2 mA 
~3-« (*-Oj=— g— 



(14) 

By combining Eqs. ([IT]), ([12)). (fT3j) and (fT3|). we eventu- 
ally arrive at the following expressions for the semiclas- 
sical diffusion coefficients, up to second order in fithuj, 

D PP (t) = 7 + ^fV(t)-7 2 ), (15) 



D qP {t) 



2j 2 mA 



(16) 



Several points are worth emphasizing. First, for undriven 
systems, the diffusion coefficients are constant in time 
and Eq. (O describes a Markovian process. This is due 
to the Ohmic regime and the semiclassical limit that we 
consider here. Incidentally, in this regime the nature of 
the initial interaction between system and bath, decou- 
pled or thermal, is not of importance. Furthermore, the 
quantum corrections D^ p and D^ p are proportional to the 
friction coefficient 7 and therefore become negligible in 
the limit of vanishing coupling. We note in addition that 
the terms 7 2 appearing in D^ }p {t) and D^ p (t), Eqs. (fT5|) 
and (fT6|) , are absent in the discussion of Ref. [12] . This 
difference can be traced back to the heuristic approach 
used in this work, which consists in taking the limit 7 — > 
for the evaluation of the diffusion coefficients. No such 
assumption is made here. 

It is instructive to check the correctness of expressions 
(fTS"]) and (|16p by computing the second moments of the 
position and momentum operators for an undriven har- 

Jq. Starting from the equation 
3), one can show that 25. 2611. 



monic oscillator, uj 2 (t) = 
for the Wigner function 

1 



(p 2 ) = 5*£ 

i»l/7 7 



D pp 




m 2 7 


m J 


m 


2m 2 A 


= 7" 


' 



1 



2A 



(3mu> 2 (3 

H - 1 2 ) , 



(17) 



(18) 



which are in agreement with the results obtained directly 
from the Langevin equation {!]). 

Quantum Smoluchowski equation. We are now in the 
position to derive the quantum Smoluchowski equation 
from the evolution equation of the Wigner function ([5]), in 
the limit of large friction. To this end, we use the method 
developed by Brinkman [27] ] and expand the Wigner func- 
tion in the basis of Weber functions D n (x) \l 



28] 



W(q,p,t) 




<Pn(q,t), (19) 



where D n (x) = 2 n / 2 e x2 l 2 H n (x) and H n (x) are the 
Hermite polynomials. The Weber functions are orthog- 
onal, f^° oo dyD n (y)D p (y) = (m\)S niP , and obey the re- 
currence relations, yD n (y) = (nD„_i(y) + D n+1 {y))/\/2 
and d y D n (y) = (nD„_i(y) - D n+1 (y))/y/2. By now 
inserting the expansion (fT9|) into Eq. ([5]) and integrat- 
ing over the momentum variable after multiplying by 
y/f3/2mexp((3p 2 /4 : m)D n (py [ / (3 /2m), one obtains the ex- 
act recurrence equation for the functions (p n (q,t), 

dt<fin + jnip n = 




d q ipn-i + (n + 1) d q ip n+ i 

1 



(3 dV 



s-¥>ti 

m oq 



1 <Pn-2 



D„ 



(20) 



In the high damping limit, or equivalently in the nonin- 
ertial limit m — > 0, terms of the form d t f n /l in Eq. (|20|) 
become negligible for n > 1 and the Brinkman hierarchy 
can be simplified to, 



dt 

<Pi 



¥>2 = 



_L dipt 

/(3m dq 
1 / dfo 




^ + (3V'(q,t)<p 



(21) 



(22) 



2jy//3m V 

+ (AV"(q,t)-mA 1 2 ) ^-^A^. 

(23) 

In the limit m — > 0, the function if2(q, t) further reduces 
to ip 2 {q,t) — (AV"{q,t) — mAj 2 ) ipo(q,t) and the set of 
three equations ([2~Tj) . ([22]) and |23|) is closed (the latter 
can be seen by multiplying Eq. ([23]) on both sides with 
\fm and by taking the m-dependence of A explicitly into 
account) . The probability density for the position of the 
Brownian particle is P(q,t) — J dpW(q,p,t) — tpo(q,t). 
Solving Eqs. (J2JJ, (JUJ) and fl23J) for tp (q,t), we then 
obtain the quantum Smoluchowski equation, 



4 



with V'(q,t) = mLd 2 (t)q and the diffusion coefficient 
D e (q,t) = (l+2AV"(q,t)). In the undriven case.EcL d24j) 
is equivalent to the equation derived in Refs. [I2|, LL4[ (in 
contrast to the one proposed in Refs. 0, [3]). Note that 
for thermodynamic reasons [13, HH, the effective diffu- 
sion coefficient D e (q,t) should be regarded as the first 
order expansion of D e (q,t) = 1/^1 — 2AV" (q,t)) (the ar- 
gument presented in Refs. 1 1 11 ] only depends on the 
diffusion coefficient and not on the explicit form of the po- 
tential). The present derivation is restricted to harmonic 
potentials and a discussion of the quantum Smoluchowski 
equation for arbitrary potentials will be given elsewhere. 

The domain of validity of the quantum Smoluchowski 
equation (|24(1 can be determined from the condition 
dt^n/l ^ [13] ■ By introducing the characteris- 
tic length scale £ 2 (t) = Dj/oj 2 (t) and using the re- 
placement d x ipx(q,t) cx ifi(£(t),t)/£(t), the requirement 
dt^ih -C <p\ leads to ui(t) -C 7, flui(t) <C kT and 
dtU)(i) <C 7 2 , which respectively correspond to high- 
friction, high-temperature and moderate-driving con- 
ditions. The constraint dt^/l ^ V2 further yields 
hj <C kT. The condition dtui(t) <C 7 2 imposes that the 
driving rate is smaller than the velocity relaxation rate 
and ensures that the nondiagonal elements of the density 
operator of the driven quantum system remain negligible 
at all times. 

Parametric harmonic oscillator. An important exam- 
ple of a driven quantum system is the parametric os- 
cillator with time-dependent frequency, ui 2 (t) — ujq + 
e 2 cos (Q,dt + <&) [H|]. Here ujq is the fixed frequency 
of the oscillator, the modulation frequency, e 2 the 
amplitude of modulation and $ an initial phase. For 
this exactly solvable system, the two functions <\>\ (t) and 
(j>2 (t) are given by Mathieu functions [15|. We men- 
tion that the properties of a classical parametric oscil- 
lator have recently been investigated experimentally in 
optically trapped water droplets [23]. For the case of 
the parametric oscillator, the moderate-driving condi- 
tion of the quantum Smoluchowski equation translates 
into efid 7 2 , showing that for fixed friction coefficient, 
the restrictions on the driving frequency become more 
stringent, the larger the modulation amplitude. 

Conclusion. We have presented a transparent and 
careful derivation of the quantum Smoluchowksi equation 
for a driven quantum system strongly coupled to a heat 
bath. Starting from the exact quantum Langevin equa- 
tion of a damped harmonic oscillator, we have combined 
a simple Green's function approach and a truncation of 
the Brinkman hierarchy in the strong friction limit to ob- 
tain the evolution equation of the semiclassical position 
distribution. Our findings confirm the results obtained in 
Refs. 12, 3] in the undriven case. On the other hand, 
we have established the range of validity of the quantum 
Smoluchowski equation for a driven system and shown 
that it restricts the driving rate to be smaller than the 
velocity relaxation rate. We have finally discussed the 



important case of the parametric harmonic oscillator. 
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